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The real-time spin dynamics and the spin noise spectra are calculated for p and n-charged quantum 
dots within an anisotropic central spin model extended by additional nuclear electric quadrupolar 
interactions (QC) and augmented by experimental data studied using identical excitation conditions. 

Using realistic estimates for the distribution of coupling constants including an anisotropy parameter, 
we show that the characteristic long time scale is of the same order for electron and hole spins 
strongly determined by the QC even though the analytical form of the spin decay differs significantly 
consistent with our measurements. The low frequency part of the electron spin noise spectrum is 
approximately 1/3 smaller than those for hole spins as a consequence of the spectral sum rule and 
the different spectral shapes. This is confirmed by our experimental spectra measured on both types 
of quantum dot ensembles in the low power limit of the probe laser. 

PACS numbers: 78.67.He, 75.75.-c, 72.25.-b 


Introduction: The promising perspective of combin¬ 
ing traditional electronics with novel spintronics devices 
lead to intensive studies of the spin dynamics of a single 
electron (n) or hole (p) confined in a semiconductor quan¬ 
tum dot (QD) n. In contrast to defects in diamonds 
d, i, such QDs may be integrated into conventional 
semiconductor devices. While the strong confinement of 
the electronic wave function in QDs reduces the inter¬ 
action with the environment and suppresses electronic 
decoherence mechanisms, it simultaneously enhances the 
hyperfine interaction between the confined electronic spin 
and the nuclear spin bath formed by the underlying lat¬ 
tice. 

Generally it is believed i i, Q, 9 that the hyperfine 
interaction dominates the spin relaxation in QDs. The 
s-wave character of the electron-wave function at the nu¬ 
clei leads to an isotropic central spin model (GSM) 9 for 
describing the electron-nuclear hyperfine coupling, while 
for p-charged QDs, the couplings to the nuclear spins can 
be mapped onto an anisotropic GSM Since the 

coupling constants for p-charged QDs are reduced com¬ 
pared to the n-charged QDs [J, [l^ , and additionally a 
large anisotropy factor A > 1 suppresses the spin decay 
of the Sz component p-charged QDs have been 

considered as prime candidates for long lived spin exci¬ 
tations in spintronics applications. 


Experimentally, however, there is evidence for com¬ 
parable spin-decay times of the Sz components 11 
in p- and n-charged QDs: hence the anisotropic GSM 
provides only an incomplete description of the relevant 
spin-relaxation processes in such systems. 


In this paper, we resolve this puzzle by investigat¬ 
ing the effect of an additional realistic nuclear electric 
quadrupolar interaction term (QG) [l^ onto the spin de¬ 
coherence. Most of the Ga and As isotopes have a nuclear 
spin / = 3/2 which is subject to a quadrupolar split¬ 


ting in electric field gradients that occur in self-assembled 
QDs by construction and couple to the quadrupole mo¬ 
ment of the nuclei 0. While previously simplified as¬ 
sumptions have been made [17Hl9l |. or the problem has 
been mapped on an effective 7=1/2 model in a ran¬ 
dom magnetic field which does not capture the full 
dynamics, we have taken into account the proper In de¬ 
pendent anisotropy and realistic strain field orientations 
estimated by a recent microscopic calculation [2lj. Al¬ 
though the short-time dynamics of p- and n-charged QDs 
are significantly different 9 , we show that the long 

time dynamics is governed by the same time scale set 
by the quadrupolar interactions in agreement with our 
experimental data presented below. 


Over the last decade, an intuitive picture for the cen¬ 
tral spin dynamics interacting isotropically with a spin 
bath via hyperfine interaction has emerged. The sep¬ 
aration of time scales [3] “ a fast electronic precession 
around an effective nuclear magnetic field, and slow nu¬ 
clear spin precessions around the fluctuating electronic 
spin - has motivated various semiclassical approxima¬ 
tions [3,0,0 1 [23-10 which describe the short-time dy¬ 
namics of the central sp in p olarization very well. As 
can be shown rigorously [23] the GSM predicts a finite 
non-decaying spin polarization [3, 10 [0 whose lower 
bound depends on the distribution function of the hy¬ 
perfine couplings and is only linked to conservation laws. 
In semi-classical theories [3, [0 it is given by a third 
of the initial spin polarization leading to a large spectral 
weight at zero-frequency in the spin-noise spectrum. The 
absence of such a zero-frequency contribution in experi¬ 
ments [nlliMIl provides strong evidence that the GSM 
is incomplete and additional interactions such as QC play 
an important role in the decoherence mechanism. 


In this work, we have employed a fullly quantum me¬ 
chanical approach, based on a Chebyshev polynomial 
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technique (GET) [32l - l34l| . to an extended anisotropic spin 
model. In order to include QC, we simulate / = 3/2 nu¬ 
clear spins. Within the GET method the largest acces¬ 
sible time scale or lowest frequency is linearly connected 
to the Ghebyshev polynomial order. All technical details 
can be found in Refs. (2^1^. 

Modelling a quantum dot: The dynamics of a single 
p- and n-charged QD is described by the Hamiltonian H 
consisting of three contributions: 


H = + Hcsm + Hqc- 


( 1 ) 


The first term represents an external magnetic field of 
strength B applied along the growth direction of the QD, 
which is defined along the z-direction. Furthermore, fj,B 
denotes Bohr’s magneton, and the occurring ^-factor de¬ 
pend on the geometry of the dots and is different for 
electrons and holes 111 . 

The coupling of the central electron or hole spin S 
to the nuclear spin bath can be casted 0 into the 
anisotropic GSM Hamiltonian Hqsm 

^ / 1 \ 
i^CSM = E + T (*^“4 + SyJl) . (2) 


Ik denotes the nuclear spin of the k-th nucleus, and N is 
the number of nuclear spins. The anisotropy parameter A 
of the spin-flip term (lOl| distinguishes between electron 
(A = 1) and hole spins, where 1 < A < oo applies de¬ 
pending on the mixture between light and heavy holes. 
Due to the enlarged Hilbert space of 2^-^+^ for / = 3/2, 
we have restricted ourselves to TV = 10 in the numerics. 
This, however, reproduces the previous results for 
N = 20 nuclear spins with I = 1/2 in the absence of the 
QG term. 

The energy scale As = J2k expected to be of 

0 ( 10 ) /xeV for electrons and approximately one order of 
magnitude smaller for holes [10| . The coupling constants 
Ak are proportional to the squared absolute value of the 
electron or hole envelope-wave function at the fc-th nu¬ 
cleus - for details concerning a realistic modelling of the 
considered set of Ak entering our numerics see Ref. (^ . 

The additional quadrupolar term [l^ in Eq. © 


N r 

Hqc=Y.<i>^ ( 4 -^ fc ) 


2 /(/+1) 


k=l 


qkv 
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(4 • nfc) - (4 • nl) 


(3) 


originates from electric field gradients in self-assembled 
QDs that couple to the nuclear electric quadrupole mo¬ 
ment and are of crucial importance for the long-time dy¬ 
namics of the central spin. The coupling constant qk is 
mainly governed by the second order derivative of the 
strain induced electric potential V [l^. The local 2:- 
direction at the A:-th nucleus is denoted by the normalized 


orientation vector which refers to the eigenvector cor¬ 
responding to the largest eigenvalue of the quadrupolar 
electric interaction tensor. The unit vectors com¬ 
plete the local orthonormal basis. 

The asymmetry parameter ry = /l A?: is com¬ 

monly neglected in the literature [17|,ll8|,l20!,l35|. A recent 
microscopic calculation of the nuclear electric quadrupo¬ 
lar couplings 2l| in self-assembled InGaAs QDs, how¬ 
ever, has found values up to ry « 0.5 depending on the In 
concentration in the QD. Therefore, we have included a 
finite 7 y = 0.5 in our calculations. 

The individual coupling constants qk are expected to 
be up to 0(l)neV [2l|, but only those qk are relevant 
for the central spin dynamics where simultaneously Ak 
is of the same order of magnitude or larger. We define 
Aq = J2k qk as a measure of relevant total quadrupolar 
coupling strength which is expected to be in the range of 
1 — lOOy^eV restricting the largest qk to gmax- The ratio 
Qr = AqjAs determines the relative strength of the QG. 

For our simulations, we generate random orientation 
vectors for each nucleus in our calculation whose de¬ 
viation angles are restricted to 9^ < 35° in accordance to 
the average deviation angle 9^ Ri 25° between the growth 
direction of the dot and the orientation vectors nj for 
Ino. 4 Gao. 6 As found by Bulutay [2l|. The coupling con¬ 
stants qk have been generated randomly from a uniform 
distribution gt/gmax S [0.5 : 1]. 

For ry = 0, Hqc partially lifts fourfold degenerate nu¬ 
clear spin states. Pinning to the growth direction, de¬ 
coherence of the central spin would be suppressed with 
increasing qk ■ A distribution of nl due to the inhomoge- 
nious strain fields [2l| favors the decoherence. Including 
a finite ry further enhances the decoherence due to the 
(5'+)^ -I- (5'“)^ term. 

The fluctuations of the transversal and longitudinal 
component of the unpolarized nuclear spin bath, re¬ 
ferred to as Overhauser field, defines the time scale 

T* = A/y/MbL+l) J2k=i ^ governing the short-time evo¬ 
lution of the central spin [v, 22| in the absence of i^Qc■ 
We have used this natural time scale to define the di¬ 
mensionless Hamiltonian H = T*H. Two factors in the 
definition of T* suggest a longer lifetime for hole spin co¬ 
herence than for electron spins: (i) the coupling constants 
Ak for holes are typically one order of magnitude smaller 
0 than for electrons, and (ii) increasing the parame¬ 
ter A > 1 to larger values suppresses flips of the central 
spin. Both factors enter the time scale linearly, yielding 
an expected lifetime increase of a factor ^ lOA for holes 
compared to electrons. However, when the spin-flip term 
in Hcsm becomes of the order of Hqc, this argument 
fails and the long time decay rate will be strongly influ¬ 
enced by the QG for p-doped QDs as we will demonstrate 
below. 

Definition of the spin-noise function: The Fourier 
transformation S(uj) of the fluctuation function 
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FIG. 1. (color online) Panel (a) shows the real time dynamics 
of a single electron confined in an InGaAs semiconductor QD 
for T* — Ins and various values of the parameter Qr = AqjAs 
without an external field. The inset shows the dependence 
of the lifetime Th, which is defined by the crossing of the 
black dashed line in panel (a) and S{t) at large times t ~ 
0(10^ — 10®) ns, on the ratio Qr on a nanosecond scale, which 
can be approximated by a power law oc ■ Panel (b) 

treats the spectral noise function S(uj) for various values of 
Qr for electron doped QDs. 

S{t) = i corresponds to 

the experimentally measured [llL l29l - l31| spectral power 
density (see below for experimental details). For very 
small probe laser intensity, all expectation values can be 
calculated using the equilibrium density operator. Hence, 
S{t) is symmetric in time, and S{u!) is given by 

/ OO poo 

S'(t)e““*dt = / S{t) cos{ujt)dt. (4) 

-OO j — OO 

From these definitions, we obtain the sum-rule 

^^(.;)=^(0) = ((^^)2)-(^Y (5) 

for the spin-noise spectrum. In the absence of an external 
magnetic field, its value is fixed to 1 /4 for a QD filled with 
a single spin. 

Since all experiments are performed in the high- 
temperature limit, the inverse temperature (3 = 0, and 
a constant density operator has been used in all numeri¬ 
cal calculations. Then the spin auto-correlation function 
S{t) also describes the spin-decay of an initially fully po¬ 
larized central spin interacting with an unpolarized 
nuclear spin bath, i. e. S{t) = {Sz{t)) /2. 

Results: For various relative QC strengths Qr, 
Fig. [1] (a) shows S{t) for electron spins (A = 1.) The 
two-stage spin dynamics is clearly visible: The initial 


FIG. 2. (color online) Spin noise of a hole doped QD in the 
(a) time and (b) frequency domain. In panel (a) S{t) is shown 
for various combinations of Qr and the anisotropy parameter 
A without an external field. Panel (b) shows the spectrum 
S{u}) for A = 4 and Qr = 1 in a varying external field B along 
the growth direction of the QD. For the B = 0 spectrum we 
supplemented a Lorentzian (dotted dashed line) with a width 
of 5 = 0.9 MHz. 


short-time decay on the scale T* to a plateau of approx¬ 
imately S'(0)/3 is only governed by the Overhauser field 
0 and not influenced by QC. Here, we have used the 
time scale of T* = 1 ns, see, for example. Ref. [l^. The 
second stage of the spin-decay is independent of the hrst 
for small values of Qr, and the decay is governed by QC. 
The shape of our curves agree remarkably with the data 
of Bechtold et. al. [l^: Qr « 0.06 — 0.1 seems to be an 
adequate choice for electrons confined in those InGaAs 
QDs. 

We have defined a second time scale Tb at which S(t) 
has dropped to the value 5(0)/6 indicated by the black 
dashed line in Fig.[T](a) (half the plateau) and have plot¬ 
ted the dependency of the lifetime Th on Qr in the inset. 
TniQr) approximately obeys a power law oc Qr 

Fig. [T] (b) shows the spin-noise spectra S{uj) for n- 
doped QDs for various Qr- The peak at around 100 MHz 
reflects the short time behavior of S{t) up to 10 ns and it 
is only slightly influenced by the variation of Qr- Since 
this peak contains approximately 2/3 of the total spec¬ 
tral weight of S{uj), the signal of the long time decay for 
electrons is expected to be a factor of 3 smaller than for 
holes. As demonstrated in Fig. [ija) the QC mainly im¬ 
pacts the low frequency peak corresponding to the long 
time decay: an increase of Qr broadens the peak width 
and induces a change of the gradient of S{uj) at interme¬ 
diate frequencies. 
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Now we focus on p-charged QDs. Since the overall QC 
strength Ag does not depend on the doping of the QD 
while Ag is decreasing by one order of magnitude when 
turning from electrons to holes, the ratio Qr is increasing 
by one order of magnitude at fixed Ag. Thus, we expect 
Th to decrease by a factor of ~ 32 when turning from 
electrons to holes. At the same time T* is increasing by 
a factor lOA, i. e. we expect the lifetime Th to be of the 
same order of magnitude for electrons and holes. 

Figure [2Ia) shows S{t) for p-charged InGaAs QDs for 
three sets of parameters A and Qr- For the conversion 
from the model parameters to the absolute time scale, we 
have assumed a reduction of Ag by a factor of 10 com¬ 
pared to the n-charged case. For fixed absolute QC pa¬ 
rameters Qk, Qr simultaneously increases also by 10, and, 
therefore, the absolute values qk are comparable to those 
used in Fig. [T] The initial decay due to the Overhauser 
field is suppressed in p-charged QDs by two effects that 
both decrease spin flips of the central spin on short time 
scales: (i) the increase of the asymmetry parameter A 
and (ii) the introduced energy splitting to the nuclei due 
to QC. Due to the lack of the short-time spin decay for 
hole spins, we define Th as S(Th) = S'(0)/e, indicated 
by the black dashed line in Fig. [2][a). For A = 4 and 
Qr = 1.0 we have determined the lifetime Th = 176 ns 
which matches the finding Th = 188 ns for electron spins 
at Qr = 0.1 extremely well. For the other parameter sets, 
the lifetimes of 400 ns {Qr = 0.8) and 740 ns {Qr = 0.6) 
are found, which are slightly larger than corresponding 
electron decay times {TH{Qr) oc Qr^^^), but still of the 
same order of magnitude. 

The spin-noise spectrum S{uj) is shown in Fig. EJb) 
for various external longitudinal magnetic field strengths 
B, X = A and Qr = 1.0 [s^. The calculated S{uj) corre¬ 
sponds to recent measurements and a nice agreement 
between our theory and the experiments is found: for 
increasing B the spectral weight, fulfilling the sum rule 
dD), is shifted from large to small frequencies. As a conse¬ 
quence the gradient oi S{u!) in the intermediate frequency 
regime uj ^ 0(0.1) MHz is increasing, which is referred 
to as a shift from an approximately Lorentzian lineshape 
for i? = 0toal// noise with increasing B as reported in 
Ref. (^ . Unfortunately, the resolution of our numerical 
investigations is limited to ~ 0.1 MHz for this parameter 
regime requiring already 6000 Chebychev polynomials. 
The linewidth of the added Lorentzian (dotted dashed 
line) at half width half maximum is 0.9 MHz, correspond¬ 
ing to the observed lifetime Th = 176 ns. Note that for 
the parameter set A = 6 and Qr = 0.8 the correspond¬ 
ing linewidth is 400 kHz, which matches the experimental 
findings of Ref. . 

For further comparison with our calculation, Fig. [3] 
shows experimentally measured spin noise spectra at a 
temperature of 5 K. The experiments were performed on 
(In,Ga)As quantum dot ensembles of similar dot density, 
in one case on average doped by a single electron per 



FIG. 3. (color online) Measured spin noise spectra at zero 
external magnetic field for electron (a) and hole spins (b) in 
ensembles of single charged QDs, measured around 890 nm 
laser excitation wavelength at 4 mW power level. The spec¬ 
tra are measured using a 650 MHz bandwidth receiver (1607- 
AC, NewFocus). Insets show measurements at 0.2 mW laser 
power using a sensitive 100MHz detector (HCA-S, Femto), 
confirming the identical zero frequency peak width. 


dot, in the other case by a single hole [ill, l3l|- The sam¬ 
ples were studied using identical excitation conditions. 
The linearly polarized light beam of a single frequency 
laser was tuned to the ground state transition energy 
maximum ll|. The laser power was reduced to 4 mW 


focused into a spot of 100 /rm diameter, giving a good sig¬ 
nal to noise ratio in 10-20 minutes of accumulation time, 
while simultaneously minimizing the laser excitation im¬ 
pact [s^ . The noise spectra are taken by a real time FFT 
using a FPGA module ll| and the spin-component is re¬ 
trieved from the noise background by interlacing the data 
at zero and 250 mT magnetic field applied in Voigt direc¬ 
tion. At 250 mT the peaked contribution to the noise due 
to spin precession is shifted out of the measured spectral 
range. 


The comparison of the electron and hole spin noise 
spectra in Figs. m a) and (b) with the calculations reveals 
that the theory qualitatively correctly predicts the shape 
and widths of the spin-noise spectra. In particular, the 
following features are worth noting: (i) The electron spin 
noise shows an additional peak around 100 MHz unveiling 
the electron’s precession in the frozen Overhauser field 
Q, as also present in Fig. dKb). (ii) Since S{u}) must 
obey the sum-rule ([5|), the low-frequency spectral weight 
of S{uj) for n-charged QD is only about 1/3 of those for 
holes. A Lorentzian fit to the low frequency components 
(/ < 35 MHz) of the experimental data confirms this 
difference in the amplitudes, (iii) a spin correlation time 
of the same order of magnitude in the long-time range 
for electrons and hole spins, as predicted by the theory. 
In the experiment this time is on the order of 400 ns, as 
estimated from the peak width at low frequencies. 


Summary: We have compared the impact of the hy- 
perfine interaction on the spin coherence in n- and p- 
charged QDs, including the nuclear quadrupolar electric 
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interaction generated by the strain fields, which provides 
an additional decoherence mechanism acting equally for 
n- and p-charged QDs. This mechanism is sufficient to 
explain the very similar long-time decay time Tjj of n- 
and p-charged QDs. On the other hand, the different 
coupling of electron and hole spins in the central spin 
part of the Hamiltonian leads to significant deviations in 
the short term dynamics, most prominently evidenced by 
the electron spin precession about the nuclear magnetic 
field. 

The samples have been provided by D. Reuter und 
A.D. Wieck, Bochum University, Germany. This work 
has been supported by the DFG and the RFBR through 
the TRR 160. 
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